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ABSTRACT 

Starting in 2012, we began an unprecedented observational program focused on the supermassive black hole 
in the center of our Galaxy, Sgr A*, utilizing the High Energy Transmission Gratings Spectrometer (HETGS) 
instrument on the Chandra X-ray Observatory. These observations will allow us to measure the quiescent X- 
ray spectra of Sgr A* for the first time at both high spatial and spectral resolution. The X-ray emission of Sgr 
A*, however, is known to flare roughly daily by factors of a few to ten times over quiescent emission levels, 
with rarer flares extending to factors of greater than 100 times quiescence. Here were report an observation 
performed on 2012 February 9 wherein we detected what is the highest peak flux and fluence flare ever observed 
from Sgr A*. The flare, which lasted for 5.6 ks and had a decidedly asymmetric profile with a faster decline 
than rise, achieved a mean absorbed 2-8 keV flux of (8.5 ±0.9) x 10~ 12 erg cm" 2 s -1 . The peak flux was 2.5 
times higher, and the total 2-10keV emission of the event was approximately 10 39 erg. Only one other flare of 
comparable magnitude, but shorter duration, has been observed in Sgr A* by XMM-Newton in 2002 October. 
We perform spectral fits of this Chandra observed flare, and compare our results to the two brightest flares ever 
observed with XMM-Newton . We find good agreement among the fitted spectral slopes (T ~ 2) and X-ray 
absorbing columns (Nh ~ 15 X 10 22 cm" 2 ) for all three of these events, resolving prior differences (which are 
most likely due to the combined effects of pileup and spectral modeling) among Chandra and XMM-Newton 
observations of Sgr A* flares. We also discuss fits to the quiescent spectra of Sgr A*. 
Subject headings: accretion, accretion disks - black hole physics - radiation mechanisms:nonthermal 



1. INTRODUCTION 

Sagittarius (Sgr) A* is the compact radio, infrared and 
X-ray source associated with the 4 x 10 6 M supermassive 
black hole at the dynamical center of our Galaxy (see, e.g., 
iMelia & Faickel 1200 It iGenzel et all I2010I) . As the nearest 
galactic nucleus, Sgr A* offers unique access to accretion 
physics on event horizon scales, and is thus a key testbed 
for theoretical modeling. On the other hand in the al- 
most 40 years since its radio identification (Balick & Brown 
119741) . a multitude of observational campaigns in the ra- 
dio/millimeter, near-infrared (NIR) and X-ray bands have es- 
tablished that Sgr A* is emitting steadily at a bolometric lu- 

1 Massachusetts Institute of Technology, Kavli Insti- 
tute for Astrophysics, Cambridge, MA 02139, USA; 
mnowak,jneilsen,fkb, houck@space.mit.edu 

2 Astronomical Institute, "Anton Pannekoek", University of Amster- 
dam, Amsterdam, The Netherlands; s.b.markoff@uva.nl 

3 Observatoire Astronomique de Strasbourg, Universite de Strasbourg, 
CNRS. UMR 7550, 11 rue de l'Universite, 67000 Strasbourg, France; 
delphine.porquet,nicolas. grosso@astro.unistra.fr 

4 School of Physics and CSPA, Monash University, VIC 3800, Aus- 
tralia; Yuri.Levin@monash.edu 

5 Leiden Observatory. Leiden University, P.O. Box 9513, NL-2300 RA 
Leiden, The Netherlands 

6 Physikalisches Institut, Unfversitat zu Koln, Zulpicher Str. 77, 50937 
Koln, Germany; eckart@phl.uni-koeln.de 

7 Department of Astrophysics, IMAPP, Radboud University Nijmegen, 
PO Box 9010, 6500 GL Nijmegen, The Netherlands; H.Falcke@astro.ru.nl 

8 ASTRON, Postbus 2, 7990 AA Dwingeloo, the Netherlands 

9 Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 
Bonn, Germany 

10 Purple Mountain Observatory, CAS, Nanjing, PR. China, 210008; 
ji@pmo.ac.cn 

11 Department of Astronomy, University of Michigan, Ann Arbor, 
Michigan, MI 48109, USA; jonmm@umich.edu 

12 Department of Astronomy, University of Massachusetts, Amherst, 
MA 01002, USA; wqd@astro.umass.edu 



minosity LboI ~ 10~ 9 LEdd, orders of magnitude lower than 
is typical for nearby lo w-luminosity active galactic nuclei 
(LLAGN; e.g.. lHo|[T999l) . Sgr A* is thus either representative 
of a distinct class of quiescent galactic nuclei lurking within 
most normal galaxies, or it is simply occup ying the extreme 
low end of the AGN continuum (see, e.g., Naga r et al.ll2005[ 
Figure 4). Placing Sgr A* into context with other objects is 
therefore an important goal, in order to correctly interpret its 
rather atypical features. 

The X-ray band is a powerful probe of the inner accretion 
flow regions of black holes. The Chandra X-ray Observatory 
soon after its launch was the first to identify Sgr A* by dis- 
cover ing a dominant, steady emission state (Baganoff et al. 
2003) which can just be spatially resolved at Chandra's sub- 
arcsecond imaging resolution. This emission can be associ- 
ated with therm al bremsstrahlu ng from near the gravitational 
capture radius (Quataert 2002; but see iSazonov et al J 120121 
for an alternative). The X-ray "quiescent state" is punctuated 
roughly daily by flares with ^ hour time scales that point to a 
source within ~ 10's of r e = GM/c 1 (the gravitational radius) 
from the black hole (Bagano ff et al. 120011: iGoldwurm et"alj 



l27il iBlianger et al.ll20T)5h iPorquet et al] 120031 120081 ). The 
flares, whose emission has been modeled with synchrotron 
or alternatively inverse Compton, are most likely caused by 
nonthermal processes (e.g., shock or m agnetic reconnection 
acceleration of electrons within a jet: iMarkoff et all 120011: 
Yuan et al j|2004t lYus ef-Zadeh et al. 2006; Dodds -Eden et all 



20091). though other me chanisms have also been suggested 



(e.g.. lLiu&Meliall 2002: Zu bovas et al.ll2012l) . 

Simultaneous monitoring campaigns with the NIR have 
established that all the X-ray flares have NIR counterparts, 
while only the brighter NI R flares (> 10 mJy ) have cor- 
responding X-ray activity dEckart et al] 12004c iGhez et al] 
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iTrap et aT]l201 ll) . The relationship of X-ray and NIR flaring 
to that in the sub-millimeter (submm) bands is still under de- 
bate, a lthough broad peaks delayed from the NIR have be en 
noted (iMarrone et al.ll200alYus"ef-Zadeh et al]l2006L 120081) . 

Many important questions persist about the nature of the 
accretion flow in Sgr A*, and particularly about the flares, 
which seem to provide a missing link with activity seen in 
other weakly accreting black holes. Specifically, in the last 
decade the Fundamental Plane (FP) of black hole accretion 
dMerloni et al.ll200l iFalcke et al.ll200l IKording et al.ll2006h 
has emerged as an important concept that links black hole 
mass and radiative output in the radio and X-ray bands for 
weakly accreting systems, with Sgr A* being an extreme low- 
luminosity example of such systems. As statistics have im- 
proved, it now appears that the X-ray flux of Sgr A* ap- 
proaches or meets the expectations from the FP only dur- 
ing its flaring state, but lies at too low an X -ray flux, rela- 
tive to its radio emission, during quiesc ence (Markoff] |2005l : 
IKording et aTj|2006t IPlotkin et alj|2012l) . As the FP radio lu- 
minosity is associated with synchrotron emission from accel- 
erated particles in compact jets, the flares may be providing 
clues about jet launching and plasma processes near the black 
hole. 

Only about two dozen X-ray flares have been detected so 
far, primarily with Chandra and XMM-Newton. Most flare 
fluxes are on the order of a factor of a few to ten times the 
quiescent flux, but a few show fluxes on the order of 100 
times greater than quiescence, and sometimes have associ- 
ated pre- or post-cursor " hiccups", i.e., weak flares close in 
time to the m ajor outburst ( Baganoff et al. 2001; Porqu et et alj 
2003, 2008). Exact flare characteristics such as spectral slope, 
which is very important for constraining models, are not 
well-determined because Chandra flares suffer pileup. Pileup 
in the Ad vanced CCD Imagin g Spectrometer-Imaging array 
(ACIS-I; iGarmire et al.l 12003) is when two or more events 
arrive in overlapping pixel regions within the same detector 
readout frame, and subsequently are either read as a single 
event with the summed energy or are discarded as a non-X- 
ray event (lDavisll200lb . Although XMM-Newton flare ob- 
servations do not suffer from pileup, its larger mirror point 
spread function (PSF) does not isolate the accretion flow as 
effectively. Prior studies had shown spectral slope differences 
between flares observed with Chandra and XM M-Newton 
dBaganoffetalJl200lt iPorquet etalj|200l 120081) . and it has 
been unclear to what extent these differences were due to in- 
strumental effects. 

As part of an unprecedented X-ray Visionary Project (XVP) 
awarded in its Cycle 13 Guest Observer Program, Chandra 
has begun the first of a total of 3 Msec of observations to be 
carried out in 2012 with the High Energy Transm ission Grat- 
ing Spectrometer (HETGS; Caniza res et al.1120051) . The main 
goal of this program is to constrain the accretion processes 
around Sgr A*, including detecting flares for the first time 
with high-resolution spectroscopy (E/AE w 200@6.4keV). 
In our first several observations, already two major flaring pe- 
riods have occurred, the first of which contains the brightest 
flare detected to date. This flare is sufficiently strong and long 
to allow us to create an individual spectrum. Here we present 
a detailed analysis of this flare and compare its spectra and 
characteristics to the two brightest X-ray flares reported in the 
last decade of Sgr A* observations. In addition we discuss the 
quiescent state, as observed in the th order of the gratings, 
and suggest a standardized method for reporting flare peaks 



and spectra to aid in characterizing the flare distribution for 
studies of the entire sample of events. 

2. OBSERVATIONS 

At the time of writing of this paper, 2012 April 1, there have 
been 39 Chandra observations of Sgr A* using its ACIS-I, to- 
taling 1.2 Msec, and now 10 using the HETG in combination 
with the ACIS-Spectroscopy array (ACIS-S), totaling 320 ks. 
The HETG is comprised of two gratings sets: the medium 
energy gratings (MEG) and the high energy gratings (HEG), 
which disperse spectra into positive and negative spectral or- 
ders. We consider only the HEG and MEG ± 1 st orders, which 
between them cover the ks 0.5-9 keV energy band. Addition- 
ally, an on-axis undispersed image is created at CCD spectral 
resolution (the th order). Compared to observations without 
insertion of the HETG, the th order efficiency ranges from 
ks 30% at 2 keV to ks 80% at 8 keV, with an average of ks 40% 
when weighted by Sgr A* 's quiescent count rate spectrum. 

The focus of our present work is a very bright flare ob- 
served with the HETGS during a 59.25 ks observation that 
began at 06:17:03 UTC on 2012 February 9 (ObsID 14392). 
The overall reduction in th order efficiency means that Sgr 
A* flares like this one are subject to significantly less photon 
pileup than ACIS-I observations conducted absent the grat- 
ings. (ACIS-S spectra also have slightly higher spectral res- 
olution than ACIS-I spectra with comparable S/N.) The dis- 
persed gratings spectra of Sgr A* flares are never subject to 
pileup. We leave the challenging analysis of the quiescent 
gratings spectra, which include significant background from 
diffuse Galactic center emission dispersed across the field of 
view, for a later date, when more data are available and we 
have a more reliable model of the background emission. 

We took two steps to isolate Sgr A*'s flare emission and 
minimize the contribution from diffuse X-ray background. 
First, we extracted spectra and lightcurves from small circular 
regions with radii of 2.5 pixels (w 1"25) centered on Sgr A*'s 
celestial coordinate^]. For the major flare discussed in this pa- 
per, we also extracted lightcurves and spectra for the ±l st or- 
der gratings (as determined by the tg_resolve_events 
tool) from long, 5 pixel wide rectangular regions centered 
on the dispersed spectra. Second, we reprocessed and ex- 
tracted CCD spectra from all existing ChandraSgr A* ob- 
servations, including ACIS-I data, to provide the best possi- 
ble characterization of the quiescent spectru m. All data were 
processed with standard CI AO v4 . 4 tools dFruscione et al.1 
2006) and calibration database v4 . 4 . 8. We selected stan- 
dard event grades (0, 2-A, 6) and applied corrections for 
Charge Transfer Inefficiency (CTI), but did not apply pixel 
randomization to the event positions. For ACIS-I and HETGS 
th order spectra, detector response matrices and effective ar- 
eas were created with the mkacisrmf and mkarf toolfl, 
with the effective areas being "aperture corrected" using the 
arfcorr tool. (This tool divides the effective area by an 
energy-dependent fraction that ranged from w 0.9-0.83 be- 
tween 2-8 keV, which accounts for the fraction of the energy- 
dependent, point source PSF within the 2.5 pixel radius source 
region.) Gratings responses were created with the mkgrmf 
(wh ich includes flux aperture c orrection) and mkgarf tools 
(see lHuenemoerder et al.ll2~01 lL for an outline of the gratings 

1 We did not re-register coordinates, since the latest Chandra data process- 
ing versions register Sgr A* to a positional accuracy of typically O'.'l. 

2 ACIS-I response matrices for ObsID 292, however, were created with the 
mkrmf tool as this observation occurred at a -1 10 C focal plane temperature. 



The Brightest Flare Seen from Sgr A* 



3 




2.8X10 4 3xl0 4 3.2xl0 4 3AxW 3.6x10" 

Time (s) 



FIG . 1 .— Top: 2-8 keV, 2-4 keV, and 4-8 keV Sgr A* lightcurves in 300 s 
bins from Chandra ObsID 14392, comprised of th order and ±l st order 
counts. Time is measured relative to the observation start: 2012 Feb. 9, 
06:17:04 UTC. Each lightcurve is fit with a constant, two Gaussian distribu- 
tions (for the precursor flares), and a Gumbel distribution. Middle: Hardness 
ratio of the 4—8 keV/2^1 keV rates, shown with the hardness ratio from the 
fits. Bottom: Close-up of the bright flare, highlighting its asymmetric profile. 

processing procedures). 

3. LIGHTCURVES 

We searched for flares in the Sgr A* lightcurves by apply- 
ing a Bayesian Blocks algorithm (Scargle 2002, priv. comm.) 
to unbinned events in the 2-9 keV band, using the imple- 
mentation from the S-lang/ ISIS Timing Analysis 
Routine j^RsiTAR). Th i s sam e method was employed pre- 
viously by Bag anoff et alj ( 12003b . We chose a detection sig- 
nificance level of 98.2%, i.e., 1 -exp(-4), for each lightcurve 
"change point". Since a flare has at least two change points, 
a rise and a decay, the overall flare significance is at least 
99.97%. Each lightcurve is then described by a series of uni- 
form rate bins (usually only one bin for the mean rate). For 
lightcurves described by multiple bins, any bin with a rate 
below the 2-cr upper bound of the lowest rate bin was consid- 

http : / / space .mit.edu/ cxc/ analysis/ SI TAR/ 



ered as "quiescent", while the remaining bins were assigned 
as "flare". Contiguous flare bins were considered to be a sin- 
gle flare, and were excised to create quiescent spectra for all 
observations. Using the Bayesian Blocks algorithm, we found 
18 flares in the ACIS-I observations and 6 flares in the HETGS 
observations, with 2 in ObsID 14392. In a future work we will 
describe the statistics and properties of the full set of flares 
detected in all Chandra observations; however, here we are 
concerned with the second flare from ObsID 14392. Both the 
mean and peak count rates from this flare were significantly 
higher than any other observed Chandra flare. 

Figure[TJpresents the full, energy -resolved X-ray lightcurve 
for ObsID 14392, comprised of th and ±l st order events 
in 300 s bins. A large flare occurs roughly halfway through 
the observation and lasts 5 ks, with two possible precursor 
flares at ~ 16,000s and ~ 22,000s. The latter flare appears 
in the Bayesian Blocks lightcurve only if we decrease the de- 
tection significance level (for detecting two change points) to 
93%. As discussed in Section |4j we create spectra for the 
mean of the large flare, but do not have sufficient statistics to 
describe the spectra at the flare peak. Determining the am- 
plitude of the flare in the lightcurve is therefore particularly 
important for assessing the flare's peak luminosity. 

The precursor flares are rather faint; however, the main flare 
is clearly visible above the background level, which is com- 
posed of the quiescent emission from Sgr A* as well as diffuse 
emission throughout the extraction regions. The main flare is 
asymmetric (Figure [JJ bottom), with a slow rise and a sharp 
decline, and can be modeled as a strong, wide (tj ~ 1400s) 
Gaussian flare followed by a weaker, narrower (er ~ 400s) 
Gaussian flare approximately 1100 s later. However, since 
our primary focus is in determining the peak brightness of the 
flare, it is useful to have this quantity as a free parameter when 
fitting the lightcurve. We accomplish this with a renormalized 
Gumbel distribution: 

/(0 = AWe ( '~' o)/T e 1 -^'-'^ , (1) 

where N pe - d k is the peak count rate, t is time, fo is the peak 
time, and r is the characteristic time scale. This provides 
a good match to the flare's asymmetry (see the red curves 
in Figure [TJ. Our full model for the lightcurve consists of 
Gaussians for the two precursor flares and a Gumbel compo- 
nent for the main flare, superimposed on a constant baseline. 
The mean flare-only count rate, calculated over the interval 
where the main flare is brighter than the background level, 
is 0.091 ±0.006 counts s" 1 ; the peak flare-only count rate is 
0.22 ± 0.2 counts s" 1 , and the characteristic time scale for the 
flare is 1 180 ± 70s. The error bars are 90% confiden ce in- 
tervals for a single parameter using the Cash statistic (Cash 
1979). Thus the flare peak rate is 2.5 ±0.3 times the mean 
flare count rate. This should be regarded as a lower limit, 
since the th order lightcurve is actually suppressed by pileup 
during the flare (see Figure|2]and the discussion below). 

There is evidence for substructure in our Bayesian blocks 
decomposition of the flare (Figure [2]): the peak of the flare 
is consistent with having a brief « 300 s dip in between two 
sharp peaks of « 100 s duration. This structure exists inde- 
pendently (albeit at lower significance) in both the th order 
and summed 1 st order lightcurves. A similar dip/short time 
scale structure was seen in the 2002 O ctober flare observed 
by XMM-Newton dPorquet et al.ll2003l) . 

We also have searched for any X-ray color differences be- 
tween the flare and non-flare intervals. First, we calculated 
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FIG. 2. — Top: For ObsID 14392, the summed th and 1 st order count rates 
(solid line) with 1<t errors (dotted line) during intervals found by the Bayesian 
Blocks algorithm (using a 98.2% significance level for each change point). 
The major flare start and stop times are 2012 Feb. 9, 14:25:50 and 15:59:51 
UTC. Bottom: For each Bayesian Block, after subtracting the mean quiescent 
level rate, the ratio of the lh order rate to the summed 1 st order rates vs. the 
summed 1 st order rates. (Error bars are Icr.) The light blue line is the expected 
correlation if the intrinsic (i.e., unpiled) th order rate is 1.7x the 1 st order 
rate, and the pileup parameter a = 1 . 

a hardness ratio (HR) as the ratio of the 4-8 keV lightcurve 
to the 2-A keV lightcurve. Zeros in the denominator were re- 
placed with the median 2—4 keV count rate. The results are 
shown in the bottom panel of FigureQ] The smooth red curve 
is the ratio of the fits to the relevant light curves, and the flare 
appears harder than the quiescent emission. A Kolmogorov- 
Smirnov (K-S) test on the extracted events indicates that the 
probability that the flare and quiescent intervals have the same 
spectrum is P = 3 x 10~ 15 (see Figure[3]for the cumulative dis- 
tribution functions, i.e., CDFs). If we consider only events 
from the ±l st order gratings spectra or the th order, the prob- 
ability that the flare and quiescent intervals have the same 
spectrum is P = 3 x 10~ 9 and P = 4.6 x 10~ 2 , respectively. We 
conclude that at the > 95% level, the flare spectrum is harder 
than the quiescent spectrum. There is no evidence for a dif- 
ference in the pre-flare and post-flare spectra, although the 
CDFs in Figure [3] do indicate that the ~ 6.6 ke V iron emis- 
sion {Baganoff et al. 2003| ISazonov et al] 120121) is relatively 
significant in the quiescent spectrum, even in this short obser- 
vation. 

To search for evidence of any color evolution during the 
flare, we examine the ratio of the th order to summed 1 st 
order count rates, as the former is most sensitive to hard X- 
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FIG. 3. — Cumulative distribution functions (CDFs) of the and ±l st or- 
der counts as a function of energy for the main flare and the quiescent interval, 
further subdivided into pre- and post- main flare intervals. The quiescent and 
flare CDFs differ at > 95% confidence. The rise in the quiescent CDF at 6.6 
keV indicates a strong contribution from the known iron emission. 

rays while the latter is most sensitive to soft X-rays. Figure|2] 
shows this rate ratio for the lightcurve sub-intervals obtained 
from the Bayesian Blocks decomposition. Because the th or- 
der is still subject to pileup, this ratio can vary even in the 
absence of color evolution. For a constant spectral shape, the 
ratio of the th order counts and ±l st order counts should scale 
as: 

(aA,)" 1 [exp(aA,)-l]exp(-A ; ) , (2) 

where A* is the incident (unpiled) counts per integration 
frame, and the odds of N piled photons being detected as a 
single "good event" is assumed to be oc cr " 1} (lDavisll2001h . 
The expected curve is shown in the bottom panel of Figure 12 
under the assumption that all piled events are recorded (i.e., 
a = 1) and that the gratings rate is «37% of the total incident 
unpiled count rate. The data are consistent with no detectable 
spectral evolution during the flare, aside from pileup effects. 
We estimate that due to the effects of pileup, our lightcurve 
measurement of the peak/mean ratio is a factor w 1 . 1 too low. 

These naive estimates, however, are subject to systematic 
uncertainties in the details of the pileup model. For the th 
order spectra, we estimate that on average, 8% of the incident 
flare photons are subject to pileup, with this fraction increas- 
ing to 16% at the flare peak. Of the detected th order flare 
events, on average 0^1% (for a = 0-1) are in fact "piled" 
events falsely registered at higher energies that harden the 
spectrum. This systematic effect is accounted for in the th 
order spectral fits below, aided by the fact that the simultane- 
ously fit HETGS spectra are not subject to pileup. Neverthe- 
less, uncertainty in pileup modeling (specifically, a) serves to 
widen the error bars on the photon index fit to the flare spectra. 

4. SPECTRA 

We next consider the spectrum for the bright flare ob- 
served in ObsID 14392, and specifically compare it to Sgr 
A* quiescent spectra. It is not completely straightforward to 
charact erize Sgr A*'s quiescent emission, as it is clearly ex - 
tended dShcherbakov & Baganoffl2010HSazonov et al.ll2012h . 
In fact, recent models of t he extended emission (i.e., 
IShcherbakov & Baganoffl20l0 ) suggest that only on the order 
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TABLE 1 

Joint Fit to Chandra Flare and Quiescent Spectra from the Sgr A* Region 
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NOTE. — The model, dust scat X TBnewX (power law+gaussian) , is applied separately to the quiescent spectrum and the flare spectrum. No gaussian 
is included during the flare. We model the th order spectrum during the flare as the sum of the quiescent emission and flare emission; for the HETG flare spectrum, 
we subtract the quiescent emission and diffuse extended emission as background. In lieu of model normalizations, we measure the integrated 2-8 and 2-10 keV 
fluxes F with the cf lux convolution model. We report both absorbed (superscript abs) and unabsorbed (superscript unabs) fluxes; the 2-10 keV luminosity L 
presumes isotropic emission at a distance of 8 kpc. F* s is the best fit model flux density at 6 keV (no errors given). Nh is the equivalent hydrogen column density, 
r is the power law index, and £i; ne , fTi; ne , and EW are the energy, 1<t width, and equivalent width of the Gaussian emission line. Errors are 90% confidence level 
for one interesting parameter. Peak flux values are derived assuming a peak/mean flux ratio of 2.5 ±0.3 (90% CL), with errors combined in quadrature. Due to 
pileup, this ratio may in fact be 10% too low (see text). 
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FIG. 4. — Left: Summed, flux corrected spectra for the inner 1 ''25 radius surrounding Sgr A*, with the fit presented in Table[T] Spectra represent: quiescent 
emission (ACIS-I, hollow blue circles; and HETG th order, solid purple squares), and the bright flare mean emission (HETG ±l st orders, solid black circles). 
Right: The same spectra (absent ACIS-I spectra) shown as detector counts s keV~' , and now including 0" 1 order spectra for the flare (hollow blue squares). 



of 1 % of the observed quiescent emission arises near the event 
horizon, as opposed to flare models where, owing to time 
scales of only thousands of seconds, almost all the emission 
is associated with the inner region. We ignore these distinc- 
tions and do not break up the quiescent emission into "point- 
like" and "extended" components, nor do we even attempt to 
"background subtract" the quiescent emission. Instead, we 
use consistent extraction regions between the quiescent and 
flare periods, and treat the flare as additional emission. Since 
our quiescent spectra include diffuse emission, we generally 
prefer to report absolute flare flux levels, rather than describ- 
ing the flare in terms of a "factor times quiescent emission." 
The mean absorbed flux density in the 2-8 keV band is the 
quantity least subject to systematic uncertainty and most use- 
ful in comparing current and prior observations. 

4.1. Methodology 

To accomplish these measurements, we created a th order 
spectrum and ±l st order spectra for the 5600s interval of the 
brightest flare in ObsID 14392. In order to isolate the actual 
flare spectrum, it is important to have a reliable characteriza- 
tion of the quiescent spectrum. For the gratings, we created a 
background spectrum by extracting a 1 st order spectrum from 
the quiescent periods of ObsID 14392. Because the th order 
spectral analysis of the flare is complicated by the presence of 



pileup, we modeled the th order flare spectrum as the sum of 
a flare component and a quiescent component. To constrain 
the quiescent component, we created quiescent spectra from 
all Chandra observations of Sgr A*. We kept the ACIS-I and 
th order ACIS-S spectra separate, but fit them with a single 
modefl This method appears justified because our K-S test 
does not indicate a difference between the pre- and post-flare 
quiescent spectra, and beca use prior studies have suggested a 
fairly stable quiescent level (Shcherbakov & Baganoff 20101) . 

We restricted our fits of the quiescent ACIS-I spectra to the 
0.5-9 keV band, the th order quiescent spectra to the 1-9 
keV band, and the gratings spectra to the 2-9 keV band. All 
combined spectra were rebinned to have a minimum signal- 
to-noise of 4.5 in each energy bin, and the gratings spectra 
were further required to have a minimum of 16 pre-binning 
channels per final energy bin. Only those bins completely 
inside the above energy ranges were included in our analysis. 

For clarity, we set out the details of our fitting process. 



4 All analyses described in this paper have been performed with 
the Interactive Sp ectral Interpretation System (ISIS; 
Houck & Denicola 2000). Spectra for each ObsID and gratings arm were 
kept separate, but like groups were combined during analysis using the 
combine_datasets function. Plots also show the combined data, and 
when showing "flux corrected" data this correction has been implemented 
using only the detector response matrices, and not the fitted models. 
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FIG. 5. — Confidence contours for Sgr A* spectral parameters at several epochs: mean unabsorbed 2-10keV flux vs. X-ray absorbing column, Nh, (left) 
and photon index vs. Nh (right). Black/blue/purple lines are for the quiescent (dotted lines) and major flare (solid lines) emission observed by Chandra in 
ObsID 14392. Brown/red /orange lines are for flare emission observed by XMM-Newton on 2007 April 4 (dashed lines) and on 2002 October 3 (solid lines) 
iPorauet et al. 2003, 2008), re-analyzed with the same spectral model as applied to the Chandra observation. Contours are 68%, 90%, and 99% significance for 
two interesting parameters. Flare flux is in units of 10~ 12 erg cm -2 s and quiescent flux is in units of 10~' 4 ergciTT 2 s -1 . 



The quiescent spectra from ACIS-I and the th order were fit 
jointly with a single model consisting of an absorbed, dust- 
scattered power-law and an iron emission line (see below for 
details). The 1 st order grating spectrum of the flare was mod- 
eled as a second absorbed, dust-scattered power-law, and the 
th order spectrum of the flare was treated as the sum of these 
two components convolved through the ISIS pileup model. 
(The pileup parameter a was left as a free parameter, but its 
error bars always spanned the full range 0-1. Again, this un- 
certainty contributes to widening the error bars on the fitted 
photon index, V.) 

Although the overall spectral model is rather simple, 
the interstellar absorption and dust scattering components 
merit further discussion. X-ray absorption is dominated 
by metals, not hydrogen, hence the fitted hydrogen col- 
umn is strongly dependent on the adopted cross sections 
and abundances (IWilms et al.l 120001) . We use the TBnew 
mo deQ develop e d from the work and abundances described 
by IWilms et al.l d2000t) . In our exper i ence, using TBnew 
with the cross sections of I Verner et al.l (11996b yields equiv- 
alent Nh values ~ 1 .5 x that o f the oft used wabs model 
( Mor rison & McCammonl 19831) . 

For all intents and purposes, given the very small Chandra 
PSF, dust scattering acts as a pure loss term, with the dust 
scattering optical depth havin g an E~ 2 dependen ce. We use 
the model dustscat (see Bagan off et al.| [2()03). which has 
an optical depth at 1 keV proportional to the X-ray absorb- 
ing column density. This proportionality has been measured 
via dust scattering hal o images and X-ray binary spectra ob- 
tained with ROSAT by lPredehl & Schmittl (119951) . who found 
r ~ 0.486 (Nh/10 22 cm" 2 ) when usi ng an analog of the wabs 
model and cross sections from iMorrison & M cCammon 
( fl983l) . Given the rough scaling between TBnew and wabs, 
we tie our dust scattering optical depth to our fitted equiv- 
alent Nh via r = 0.324(N H /10 22 cm" 2 ). The implied extinc- 
tion (from the correlations of Predehl & Schmitt 1995) is then 
Av ~ N H Bnew /2.69 x 10 21 cm" 2 . However, these dependencies 
have not been revisit ed with modern absorption or dust scat- 
tering models (e.g., Xiang et al. 201 1, and references therein) 
using instruments capable of both imaging halos and making 

5 http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/index.html 



direct measurement of metal absorption edges (i.e., Chandra- 
HETG and XMM-Newton -RGS), and must be treated as hav- 
ing a certain degree of systematic uncertainty. However, in 
regard to the values presented in Table Q] adopting a ratio of 
r to N H /10 22 cm" 2 that lies between 0.243-0.486 alters the 
implied X-ray fluxes by only ±5%, our fitted Nh values by 
±10 22 cm~ 2 , and our fitted photon indices by ±0.05. 

4.2. Results 

Our fit results are presented in Tableland are shown in Fig- 
ure!?] Consistent with the results of the K-S test for the single 
ObsID 14392, the summed quiescent spectrum is significantly 
softer than the mean flare spectrum, with no overlap between 
the 90% confidence level error bars for their photon indices, 
r. The iron emission line at 6.63 keV in the quiescent spec- 
trum has an equivalent width (EW) of «780 eV. If, instead of 
a power-law, we fit the quiescent continuum with a thermal 
plasma model that already includes ionized iron emission at 
a slightly higher energy than above (vmekal), we then find 
evidence for a 120 eV EW line at 6.3±[} 2 keV, consistent with 
Fe Ka fluorescence, and consistent with the previous sugges- 
tion of such a line by Sazonov et al] (120121) . This model pro- 
vides a similarly good fit (x 2 /DoF= 264/256) with plasma 
temperature kT = 2.7^q | keV and a slightly smaller X-ray ab- 
sorbing column density Nh = (1 1 -8^o g) x 10 22 cm" 2 . There is 
good agree ment with previous studies of the quiesce nt spec- 
trum (e.g., Bag anoff et al.ll2003l:ISazonov et alfcoTH) . despite 
the fact that these authors used larger extraction regions (1" 
.5 radius) and very different background subtraction methods. 
We note that there is excellent agreement between the th or- 
der quiescent spectrum and the ACIS-I spectrum. 

By any measure, the bright flare mean emission observed 
in ObsID 14392 is significantly brighter than the quiescent 
emission, and assuming that the peak spectrum has the same 
spectral shape as the mean spectrum, the peak emission is 
a factor of 2.5 times brighter still. The mean, absorbed 2- 
8 keV flux is higher than any Sgr A* mean flare flux observed 
with Chandra, and almost identical to the mean flux of the 
brightest S gr A* flare ever ob s erved by XMM-Newton (2002 
October; iPorquet et ail I2003L 120081 and Section |5]). With 
a flare duration of « 5.6ks (compared to » 2.8 ks for the 
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2002 October flare), the flare's absorbed 2-8 keV fluence is 
(4.7 ± 0.5) x 10" 8 erg cm" 2 , and its emitted intrinsic energy in 
the 2-10keV band is approximately 10 39 erg. As is evident in 
Table [T] the further we extrapolate beyond the well-measured 
2-8 keV band, and if we consider unabsorbed instead of ab- 
sorbed fluxes, the greater the uncertainty becomes both statis- 
tically and systematically. 

Whereas the mean, absorbed flux is extremely well- 
constrained, the unabsorbed fluxes have strong dependencies 
upon the fitted X-ray absorbing column. Confidence contours 
for quiescent and mean unabsorbed flare flux vs. equivalent 
neutral column are presented in Figure [5] This figure also 
shows confidence contours for fitted photon indices, T, vs. 
equivalent neutral column. Unsurprisingly, there are strong 
correlations between indices and columns, with harder pho- 
ton indices being associated with lower columns. There are, 
however, two important points to note: there is a good consis- 
tency between the fitted columns for the quiescent and flare 
spectra (with the latter allowing a wider range of values ow- 
ing to the poorer statistics), and there is a very clear separation 
between the index/column contours for the quiescent and flare 
spectra. The flare spectrum is significantly harder; however, 
for most values of X-ray absorbing column the flare photon 
index is somewhat softer than implied by previous Chandra 
measurements of a bright flare that indi cated T = 1 ±0.8 (90% 
confidence level; Baganoff et al. 2001). Absent the use of the 
HETGS, however, there is a question as to what extent this 
prior result was influenced by pileup. This prior result also 
used an earlier version of the Chandra calibration and did not 
include the E~ 2 dependence of dust scattering losses to the 
spectrum. 

5. COMPARISON TO FLARES OBSERVED BY XMM-Newton 

Our best-fit mean flare photon index of T = 2^ 6 is in good 
agreement with previous results obtained for the brightest 
(2002 October 3) and second brightest (2 007 April 4) Sgr A* 
flares ever observed by XMM-Newton ( IPorquet et alj|2003l 
2008), from 22 observations centered on Sgr A* which totaled 
1.1 Msec through 2009 April. In order to carefully assess 
the degree to which the properties of the HETGS observed 
flare are comparable to the brightest XMM-Newton observed 
flares, we have re-fit the 2-10 keV spectra of the 2002 Octo- 
ber 3 and 2007 Apr il 4 flares — using the same data files as 
IPorquet et alj (|2008) — with the identical absorbed/scattered 
power-law modeQ presented in Section|4] 

For the 2002 October 3 flare we find a photon index T = 
2.3 ± 0.3, equivalent neutral column of Nh = (I6.I+22) x 
10 22 cm" 2 , and an unabsorbed 2-10keV flux of (26.0^'j) x 
10~ 12 erg cm" 2 s" 1 . This flare's mean absorbed 2-8 keV flux 
is (7.7 ±0.3) x 10~ 12 ergcm" 2 s" 1 , which for its 2.8 ks dura- 
tion corresponds to an absorbed 2-8 keV fluence of (2.2 ± 
0.1) x 10~ 8 erg cm" 2 and an intrinsic emitted energy in the 2- 
lOkeVband of5.3 x 10 38 erg. The phot on indices found here 
are very similar to those reported in lPorquet et alj d2008l) with 

6 IPorquet etal] 120031) analyzed the brightest XMM-Newton Sgr A* flare 
with a dust scattering model that presumed a fixed Ay = 30 (i.e., the scatter- 
ing optical depth was not tied to the fitted N h) and used the tbabs absorp- 
tio n model an d abundances of Wilms et al. (2000) with the cross sections 
of IVerne r et al. ( 1996). Porauet et al. (2008), in order to compare to other 
works, instead analyzed the two brightest XMM-Newton flares with a dust 
scattering model wi th a fixed A v = 25, and used wabs w ith cross sections of 
Morrison & McCammon (1983) and the abundances of Anders & Grevesse 



a AT of only +0. 1 . The Nh values differ due to the different 
abundances and cross-sections assumed in this work. 

We find a photon index of T = 2.4^ j for the 2007 April 
4 flare, an equivalent neutral column of Nh = (16.3^ 6) x 
10 22 cm" 2 , and an unabsorbed 2-10keV flux of (16.8!^) x 
10~ 12 erg cm" 2 s . The mean absorbed 2-8 keV flux is 
(4.8!„ 2 ) x 10" 12 erg cm" 2 s" 1 , which for the flare's 2.9 ks du- 
ration corresponds to an absorbed 2-8 keV fluence of (1.4 ± 
0.1) x 10~ 8 erg cm" 2 and an intrinsic emitted energy in the 2- 
lOkeV band of 3.5 x 10 38 erg. 

Figure|5]also presents confidence contours of X-ray absorb- 
ing column vs. 2-10 keV unabsorbed flux and photon index, 
r, for these two XMM-Newton observed flares. We see that 
the 2002 October and 2012 February events appear to be "twin 
flares" in all respects, and that aside from having a lower flux, 
the 2007 April flare otherwise appears identical to these two 
extremely bright flares. The XMM-Newton observed flares 
have slightly larger best-fit values for the equivalent neutral 
column; however, there is a high-degree of overlap among the 
error contours. Both XMM-Newton observed flares are also 
slightly shorter in duration, lasting « 3 ks, with the brightest 
flare lightcurve also showing a brief dip near its peak, de tected 
in all three EPIC instruments (IPorquet et al.ll2003Ll2008h . 

6. DISCUSSION 

Our Chandra -HETGS observation of Sgr A* taken on 2012 
February 9 (ObsID 14392) exhibits a flare with the highest 
peak flux and fluence seen from this source. Remarkably, 
it is bright enough to allow the extraction of a pure flare 
gratings spectrum. Our comparative analysis indicates that 
in many ways, the bright Sgr A* flares observed by XMM- 
Newton in 2002 October and 2007 April are spectral twins 
to this flare. Depending upon how one defines and measures 
Sgr A* 's quiescent flux and in how one extrapolates the mean 
flare flux to a peak flux, for both the 2012 February and 2002 
October flares, the ratio of peak flux to quiescent flux is at 
least a factor of order 130. If o ne adopts the suggestion of 
lShcherbakov&" Baganoff (120101) that only « 1% of the ob- 
served quiescent flux is from the central regions (where the 
flare likely originates), then this factor is more plausibly of 
order 10 4 ! 

Given systematic uncertainties in extrapolating unabsorbed 
fluxes and defining the quiescent flux associated solely with 
the Sgr A* point source, however, we suggest that a less am- 
biguous set of reported values are the flare's mean absorbed 
2-8 keV flux, its absorbed 2-8 keV fluence, and the ratio of its 
peak rate to mean rate, each measured as values above the qui- 
escent level and aperture corrected for the instrument's PSF 
So long as instrumental spectral extraction regions are consis- 
tent for quiescent and flare spectra, these values will be well- 
defined. (However, due to the short time scale sub-structure 
in the flare, as see n in Figure [2] and previously reported by 
IPorquet et alJ2003l the "peak" flux value may actually be dif- 
ficult to define precisely.) The 2-8 keV band is above the 
range of the most severe X-ray absorption in Sgr A* and is 
well-covered by the three soft X-ray instruments best capable 
of observing Sgr A*: Chandra, XMM-Newton, and Swift. 

The bright flare is asymmetric, with a slower rise than 
decay. This may be due to unresolved sub-structure: the 
Bayesian Blocks decomposition shows evidence of complex 
structure near the flare peak, and a fit with two Gaussian pro- 
files, with different widths and offset from one another, works 
well. (The fitted Gumbell profile is a convenience that allows 
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us to easily calculate the flare peak/mean ratio.) Alternatively, 
if this is the flare's intrinsic profile it is quite different than the 
"fast rise, exponential decay" of many different types of tran- 
sient phenomena. 

The total emitted 2-10keV energy of the flare, which is 
of 0(lO 39 erg), requires conversion of at least 10 19 g of rest 
mass into energy, presuming a 10% conversion efficiency. 
There is no agreed upon mechanism for flare energization, 
with suggestions having ranged from magnetic reconnec- 
tion dMarkoff et alJl200UlYuan et al.ll2003l:lDodds-Eden et al l 
20101) to tidal d i sruption of asteroid s (ICadez et alJ 120081: 
Kostic et al.l2009tlZubovas et al.l2012l) . Regardless of the en- 
ergization mechanism, with the results from the prior XMM- 
Newton observations it appears that at least these extremely 
bright flares require an emission mechanism that produces a 
moderate (i.e., not very hard) photon index, T « 2. 

As we have shown, much of the uncertainty on the spectral 
slope is systematic, and depends on assumptions made about 
absorption and scattering. Reasonable assumptions about the 
scaling between dust scattering and absorption have only a 
small systematic effect (AT ~ 0.05), but the dependence upon 
fitted Nh is more pronounced (i.e., Figure[5]). Although there 
is a good consistency between the equivalent neutral column 
for the quiescent and flare spectra, it is possible to alter the fit- 
ted Nh by assuming different metal abundances, or metal de- 
pletions in dust grains, or different cross sections , etc. Based 
on the correlation from iPredehl & Schmittld 19951) and our as- 
sumptions about A^h and Ay outlined in Section [4] the min- 
imum reasonable X-ray absorbing column for the quiescent 
spectrum implies an extinction of Ay ~ 40. This is somewhat 
at odds with prior estimates of Ay: a discussion of the Galac- 
tic cejite£_oj)tical_and infrared extinction curves can be found 
in iFritz et all (1201 lb . who suggest a value of Ay ss 33 (in 
pa rt based upon the X -ray observations of Sgr A* dis cussed 
in lPorquet et al.ll2~003l) . The relationships derived by Watson 
J20T1 and lGuver & Ozel J2009h . A v « N H /2.2 x 10 21 cirT 2 , 
imply a similar value for the extinction towards Sgr A* . How- 
ever, the methods used to derive jVh and Ay in these works are 
very heterogeneous, and it is unclear if their scalings can be 
extrapolated to Ay ~ 40. It is therefore imperative that scaling 
between A^h, extinction, and dust scattering optical depth be 
revisited with modern models using consistent cross sections 
and abundances, and with modern high spectral and imaging 
resolution observations. For the moment, we conclude that 
the extinction, dust scattering optical depth, A^h, and (by ex- 
tension) the spectral properties of Sgr A* are still subject to a 
certain degree of systematic uncertainty. 

Despite these uncertainties, we find similar spectral proper- 
ties for the brightest Chandra and XMM-Newton flares (T ~2 



and A^h ~ 15 X 10 22 cm" 2 ). Weaker flares could in princi- 
ple have harder spectra than bright flares, but their spectral 
properties ar e not yet yet strongl y constrained (as shown and 
discussed in Porquet et alJ 120081) . However, there is signif- 
icant cause for optimism, as the observations discussed in 
this work represent slightly more than 10% of the Chandra- 
HETGS observations of Sgr A* that will occur in 2012. We 
anticipate that this program will detect over three dozen flares, 
with perhaps one or two more with amplitudes comparable to 
or greater than that of the 2012 February flare. As we have 
demonstrated that we can extract Chandra gratings spectra for 
individual bright flares, the XVP program should provide un- 
precedented constraints on the spectral properties of faint and 
moderate flares, which will allow us to determine how the 
physics of flares scales with their luminosity. 

For the first time, the spatially resolved Chandra studies 
will produce Sgr A* flare spectra that are either absent of de- 
tector pileup or that have pileup strongly constrained (as is the 
case here) owing to the simultaneous 1 st order gratings spec- 
tra. Our understanding of the physics of flares will be greatly 
enhanced by the fact that many of the upcoming observations 
will be performed with simultaneous multi-wavelength obser- 
vations. The Chandra XVP program will offer us the unique 
opportunity to study the physics underlying accretion onto Sgr 
A* and other quiet galaxies. 
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